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O ■ Abstract 
(S, 

We study an extended system that without noise shows a spatially homogeneous state, but when 
submitted to an adequate multiplicative noise, some noise-induced patterns arise. The stochastic 
resonance between these structures is investigated theoretically in terms of a two-state approx- 
jy^ , imation. The knowledge of the exact nonequilibrium potential allows us to obtain the output 

signal-to-noise ratio. Its maximum is predicted in the symmetric case for which both stable attrac- 

. - 

tors have the same nonequilibrium potential value. 
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I. INTRODUCTION 



During the last few decades a wealth of research results on fluctuations or noise have 
lead us to the recognition that in many situations noise can actually play a constructive 
role that induces new ordering phenomena. Some examples are stochastic resonance in 



zero-dimensional and extended systems 
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a, 



noise-induced transitions jjjj, noise- 
noise- sustained 



induced phase transitions [8J, |9|, noise- induced transport 
patterns 14 . Il5 , Ir]] . noise-induced limit cvcle 17 1, etc. 

The phenomenon of stochastic resonance (SR) — namely, the enhancement of the output 
signal-to-noise ratio (SNR) caused by injection of an optimal amount of noise into a nonlin- 
ear system — stands as a puzzling and promising cooperative effect arising from the interplay 
between deterministic and random dynamics in a nonlinear system. The broad range of 
phenomena — drawn from almost every field in scientific endeavor — for which this mech- 
anism can offer an explanation has been put in evidence by many reviews and conference 
proceedings. See Ref. [l| and references there to scan the state of the art. 

Most of the phenomena that could possibly be explained by SR occur in extended systems: 
for example, diverse experiments were carried out to explore the role of SR in sensory 
and other biological functions [18] or in chemical systems These were, together with 
the possible technological applications, the motivation to many recent studies showing the 
possibility of achieving an enhancement of the system response by means of the coupling of 
several units in what conforms an extended medium P,!^ 



j, lij, L«a, iu, 



l20| , or analyzing the possibility 




of making the system response less dependent on a fine tuning of the noise intensity, as well 
as different ways to control the 

In some previous papers 
in extended systems for the transition 
concept of nonequilibrium potential (23, 



renomenon |21|, |2J 

we have studied the stochastic resonant phenomenon 
aetween two different patterns, and exploiting the 



24j. The nonequilibrium potential is a special Lya- 



punov functional of the associated deterministic system which for nonequilibrium systems 
plays a role similar to that played by a thermodynamic potential in equilibrium thermo- 
dynamics 2^. Such a nonequilibrium potential, closely related to the solution of the time 
independent Fokker-Planck equation of the system, characterizes the global properties of 
the dynamics: that is attractors, relative (or nonlinear) stability of these attractors, height 
of the barriers separating attraction basins, and in addition it allows us to evaluate the 
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transition rates among the different attractors. 

In this work we analyze a new aspect of such a problem studying SR between two noise- 
induced-patterns jlJ . that is: the same noise source that induces and supports the 
existence of the patterns is the one that induces the transitions among them, and produces 
the stochastic resonant phenomenon. Some related work correspond to the misleadingly 
called double stochastic resonance |25], as well as to another previous work 2o| related with 

nir 

noise-induced phase transitions |8|, |9( . In both cases the authors have only resorted to a stan- 
dard mean-field approach, while here we obtain the exact form of the noise-induced patterns 
(stable and unstable ones) as well as the complete form of the nonequilibrium potential. 
In this way we can obtain the transition rates and clearly quantify the SR phenomenon by 
means of the SNR. 

The organization of the paper is as follows. In section |H] we present the model and 
formalism to be used. After that, we discuss in section ITTT1 the stochastic resonance between 
the noise-sustained structures. Finally, we present in section H3 some conclusions and future 
perspectives. 

II. THE MODEL 

We consider a one-dimensional system described by the following deterministic equation 

d t <f ) (x,t) = d x (D(<f ) )d x <P) + F(<p), (1) 
that can be written in a variational form as 



where the potential 



L/2 ^ -I 



V[<P}= dx{- d<f>'D(4> l )F(<f>') + -[D(4>)d x 4>V\ (3) 

is a Lyapunov functional (while D{4>) > 0) for the deterministic dynamics and it is essentially 
the logarithm of the probability density of configuration when Eq.(£Q) is perturbed by an 
additive source of spatiotemporal white noise. 

The starting point of our stochastic analysis will be Eq. $!]) with an additional multi- 
plicative noise, in the Stratonovich interpretation, given by 

d ^ = -mW) +g(mxA (4) 



where £ is a Gaussian noise with zero mean and correlation (£(x, t)^(x', t')) = 2e 2 5(x — 
x')5(t — t'), being e 2 the noise intensity. For the coefficient of the noise term, g(<fi), we adopt 



in order to guarantee that the fluctuation-dissipation relation is fulfilled 

As we are considering the Stratonovich interpretation, the stationary solution of the 
associated Fokker-Planck equation can be written as j3] 

P st [0] ~ exp(-K///e 2 ), (6) 
where the effective potential V^//[0] is given by 

rL/2 

Veff[4>] = V[4>] - A / dxhxD{4>). (7) 

J -L/2 

Here A is a renormalized parameter related to e through A = e 2 / (2Ax) in a lattice discretiza- 
tion, where Ax is the lattice parameter 

The extremes of V e ff correspond to the (stationary) noise-sustained structures (f> s t. They 
can be computed from the first variation of V e ff(4>) respect to <fi equal to zero, that is 

rL/2 

5V eff [<f> st \ = - / D{4>)[d x {D{4>)d x <t>) + F eff {<t>)\5<t>{x) dx = 0, (8) 

J-L/2 <t>=<i>st 

where 

is the effective nonlinearity which drives the dynamics of the noise sustained patterns. 

We will consider a finite system, i.e. limited to the region —L/2 < x < L/2; and assume 
Dirichlet boundary conditions, that is <p(±L/2) = 0. In addition, we consider the case of a 
monostable dynamics in absence of noise 

F = —(f) 3 + b(j) 2 , (10) 

and we adopt a model of field-dependent diffusivity which induces an effective bistable 
dynamics. In particular we have chosen 
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(D , h > 0), that corresponds to have a larger diffusivity in low density (low 0) regions and 
a lower diffusivity in high density (large 0) ones. With this functional form, F e ff(<j)) in Eq. 
© results 

F eff = -0 3 + 60 2 - ^ = 0(0 - 0x)(0 2 - 0), (12) 

where 0i 2 depend on parameters, in particular on the control parameter A. It is worth 
noting here that in the deterministic problem (A = 0) the reaction term is monostable 
while, as we increase the noise intensity, the effective nonlinear term F e ff becomes bistable 
(within the interval < A < b 2 Do/(8h)) and finally, for A > b 2 Do/(8h) becomes again 
monostable (reentrance effect). Our choice of F and D is one among a plenty of different 
forms for the diffusivity leading to a transition from monostable to bistable and inducing 

28\ . that corresponds exactly to the 



the SR phenomenon (see for instance the one used in 
inverse of the present diffusion coefficient, i.e. -D(0) = D (l + hcf 2 )). Density- dependent 
diffusivities arise in a large variety of systems modeled by reaction-diffusion equations 
In biology, for instance, population dynamics is usually driven by a diffusivity that depends 
on the local population j30|. We can also find examples in physics, a couple of them are in 
polymer physics (where the diffusion can abruptly drop several orders of magnitude at the 
gelation point js^l) and in diffusion of hydrogen in metals j^. 

A remarkable point is that = is always a root of F e ff = (see Fig. 1). This implies 
(from Eq. (|5|)) that <f)(x) = is an extremum of V e ff[<f)] for all values of A. In what follows 
we will call this pattern 0o- 

In order to obtain the non uniform extremes of the potential (and also of the density 
probability) we must (numerically) solve 

±(D(<p st )4-<l>st) +F e// (0 st ) = O, (13) 



dx V dx 

for the stationary regimen profiles <p st (x). This approach allows us to found both, the stable 
and unstable solutions. To analyze their stability we need to calculate 5 2 V e ff, that defines 
a Sturm-Liouville problem, with orthogonality weight D{<p st ). From that analysis it results 
that 0o (defined before) is stable for A > 0, and in the bistability region we have two 
nonhomogeneous symmetric patterns: one unstable <p u (saddle) and one stable S . The 
typical form of these patterns is illustrated in Fig. 2. 

In Fig. 3 we show V r e //[0st] vs. A, evaluated on the different stationary patterns. We define 
A c as the value of A at which we have symmetrical stability, i. e. where K//[0o] = K//[0s]- 




FIG. 1: Form of the nonlinearities for the deterministic case (A = 0), bistable case (A = 0.8) 
and a monostable case (A = 1.2) in the reentrance region. The vertical scale was changed in the 
deterministic case in order to clarify the figure. The parameters used are: Dq = 1, h = 1/2 and 
6 = 2. Note that = remains as a root in all cases. 

From this figure it is apparent that, similarly to previous cases j^j], for increasing values of 
A both nonhomogeneous solutions (the stable and the unstable one) coalesce and disappear 
in agreement with the above argument about the reentrance effect. 



III. STOCHASTIC RESONANCE BETWEEN NOISE-SUSTAINED STRUC- 
TURES 

We are interested in the stochastic resonance phenomena occurring in the above described 
system. For a window of noise intensity the effective dynamics of the system is bistable, 
corresponding to a noise-induced nontrivial dynamics. We will resort to the so-called two- 
state approximation j^, all details about the procedure and the evaluation of the SNR 
could be found in j^. We consider now that the system is subject, in the adiabatic limit, to 
a time periodic signal of the form b = b + S(t) where S(t) = Ab sin(u;o£). 

Up to first-order in the amplitude Ab (assumed to be small in order to have a sub- 
threshold periodic input) the transition rates Wi take the form 

W\{t) = — ai Ab sin(u;o t), 
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FIG. 2: Noise-sustained stationary patterns of the problem. We show <po = 0, the stable homoge- 
neous pattern; and both nonhomogeneous patterns: the unstable (saddle) <p u and the stable one 
4> s . Here we have A = A c ~ 0.8, while the other parameters values are the same as in Fig. 1 



W 2 {t) = jj, 2 + a 2 Ab sm(u Q t), 



(14) 



where the constants [i\^ 2 and a\ )2 are obtained from the Kramers-like formula for the tran- 
sition rate j^] 

' det VefM 1 1/2 
J det V eff [(p u }\_ 

Here A + is the unstable eigenvalue of the deterministic flux at the relevant saddle point ((f) u ) 
and 



exp[-(V; // [0 u ]-^ // [0 l ])/ e 2 ]. 



(15) 



^1,2 
«1,2 



W lj2 \s(t)=a 
dW h2 , 

S(t)=0- 



dS(t) 



(16) 



These results allows us to calculate the autocorrelation function, the power spectrum and 
finally the SNR. The details of the calculation were shown in Ref. jjj. For the SNR, and up 
to the relevant (second) order in the signal amplitude Ab, we obtain 



R 



4/Xi/i 2 /il+yU2 



4e 2 ^ + \i 2 



where 



L/2 r-4>s{x) 

dx / D( 

-L/2 J<j> 



k'2 



(17) 



(18) 




FIG. 3: Nonequilibrium potential V e ff[<p s t], as a function of A, evaluated on the stationary patterns: 
curves correspond to stable (<fi s ), homogeneous (<fo) and unstable (c/> u ) patterns. The arrow indicate 
the point where K//[0o] = ^e//[^s]j corresponding to A = A c « 0.8. 

gives a measure of the spatial coupling strength. In our case (fi = and 

$ = D J^ dx j^M _ arctan[y^^(x)]| (ig) 

In Fig. 4 we show the SNR as a function of parameter A (which is proportional to e 2 ). The 
existence of the typical maximum is the characteristic fingerprint of SR. For a window of 
noise intensity values, the system enhances the output to the input periodic signal. We see 
that the maximum SNR occurs at the symmetric situation, that is at A = A c . 

A similar behavior is observed in general for a wide range of values for h and Dq com- 
patible with a bistable effective dynamics. In particular, A c is a monotonically decreasing 
function of h, as we show in Fig. 5. For a given value of h, a numerical analysis of Eq. (fTTj) 
indicates that the maximum of SNR take place at X c (h). Note that, for a given value of A, 
h appears as a additional control parameter that allows a fine tuning of the symmetrical 
condition. Finally, in Fig. 6 we show R c = R(X C ) vs. h in the range of values where Kramer's 
formulae applies |3fij . 
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FIG. 4: Signal-to-noise ratio vs. A. Here 6q = 2, while other parameters remain unchanged. 
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FIG. 5: A c vs. h parameter of diffusivity. For small h values A c , and hence the noise intensity, 
increase monotonically. 

IV. CONCLUSIONS 



The study of SR in extended or coupled systems, motivated by both, some experimental 
results and the technological interest, has recently attracted considerable attention [2, La, m 
BBQ. m SOme previous pape, flQQQ we have studl ed th e SR pae_L the 
;ransition between two different patterns, exploiting the concept of nonequilibrium potential 
2^1 Q| • In this work we have analyzed the SR phenomenon in an extended system from a 
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FIG. 6: Signal-to-noise ratio at A c vs. h. We can see a saturation phenomena as h decreases. 



different point of view, that is studying SR between two noise-induced patterns [14, 1! 
Some related work correspond to the misleadingly called double stochastic resonance 



as well as to a previous work 



[25]. 



that is tightly related to noise-induced phase transitions 



0. In both cases the authors have resorted to a standard mean-field approach. Here 
we adopt a different approach, obtaining numerically the exact form of the noise-induced 
patterns (both the stable and unstable ones) as well as the analytical expression of the 
nonequilibrium potential. In this way we were able to obtain the transition rates and clearly 
quantify the SR phenomenon by means of the SNR. 

We have seen that the a nonhomogeneous spatial coupling, through density-dependent 
diffusivity, changes the effective dynamics of the system and, in agreement with 36[, that 
such nonhomogeneity could contribute to enhance the SR phenomenon. The form of the 
patterns, position of the attractors, barrier's high, explicitly depend on the noise intensity. 
We have found that there are ranges or windows of noise intensities where the phenomenon 
could arise (reentrance). 

By considering the adiabatic limit and exploiting the two-state approximation we have 
theoretically predicted the occurrence of SR between those noise-sustained patterns. It is 
worth here remarking that it is the same noise source the one that sustains the patterns 
and induces SR for transitions among them. The maximum of the SR response occurs in 
the symmetric case, in agreement with the results found in j^, The SR phenomenon 
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is robust respect to variations of the h parameter of diffusivity, and when h decreases the 
SNR maximum increases and shifts toward higher A values. The last fact follows from the 
associated shift of the noise-induced transition to larger noise intensities which take place 
in the spatially uncoupled associated system (i.e. the 0-d system resulting from suppressing 
the gradient term in Eq.(jHJ). 

The consideration of more general forms of couplings in many component systems will 
allow us to analyze SR between noise-induced patterns in activator-inhibitor-like systems. 
We will also study, within the present framework, the competence between local and non- 

nn 

local spatial couplings |4|,|f|, etc. These aspects, together with Monte Carlo simulations of 
the different cases, will be the subject of further work. 
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